Entangling homogeneously broadened matter qubits in the weak-coupling cavity-QED regime 
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In distributed quantum information processing, flying photons entangle matter qubits confined in cavities. 
However, when a matter qubit is homogeneously broadened, the strong-coupling regime of cavity QED is typ- 
ically required, which is hard to realize in actual experimental setups. Here, we show that a high-fidelity en- 
tanglement operation is possible even in the weak-coupling regime in which dampings (dephasing, spontaneous 
emission, and cavity leakage) overwhelm the coherent coupling between a qubit and the cavity. Our proposal 
" " -." , enables distributed quantum information processing to be performed using much less demanding technology 

j ■ than previously. 

^ ! PACS numbers: 03.67.Bg, 42.50.Ex, 03.67.Lx 

Distributed architecture is a promising approach for realizing scalable quantum computation |H-@]. Elementary nodes com- 
posed of a few qubits are networked to achieve scalable quantum computation. The node separation can potentially suppress 
decoherence induced by uncontrollable interactions between qubits. Moreover, since the nodes are spatially separated, individual 
7—1 \ qubits can be easily addressed by the optical field. 

A critical operation for realizing distributed quantum computation is the entanglement operation (EO) lU-IH]. To construct 
C"| | an entire network, qubits in distant nodes have to be coupled by EOs. Most EOs are based on photon interference, and the 
successful execution of an EO can be heralded by detecting a photon at the target port. This approach has been experimentally 
+i> ' demonstrated using an ion trap system ifioll . If the EO fails, the two qubits involved should be initialized, which risks destroying 
the entanglement of other qubits generated by previous EOs. Although EOs typically have such probabilistic properties, previous 
studies have revealed that only polynomial steps are required to construct large entangled states B fTTUTill . such as a cluster state 
qh! lfl5tl . Moreover, by introducing a quantum memory to each node, EOs can be repeatedly performed until they are successful 

■ without destroying prior entanglement iflrl . 

EOs involve optical excitations of matter qubits. However, the excited states are inherently noisy and significantly degrade 
the target entanglement. For example, nitrogen vacancy (NV) centers in diamond have promising properties such as a long 
coherence time at room temperature and optical addressability. Entanglement between an NV center and an emitted photon 
has been demonstrated at a low temperature of about 7 K lfl7ll . However, at room temperature, this otherwise attractive system 
suffers from strong environmental dephasing originating from interactions with phonons when the system is optically excited. 
Consequently, it acquires a large homogeneous broadening of the order of THz 01 811 . Therefore, in such an approach, NV centers 

■ can be used for distributed quantum computation only at low temperatures. 

| One way to overcome homogeneous broadening is to employ high-Q cavities. Previous theoretical proposals of EOs require 
. strong coupling between a matter qubit and the cavity when the matter qubit has large homogeneous broadening IU9l - [2lll . 
t— I ' However, despite rapid advances in cavity fabrication technology, it is still very difficult to experimentally generate strong 
. coupling between a matter qubit and a high-Q cavity. To realize distributed quantum computation, it is thus essential to examine 
• i-h ' the possibility of performing an EO in the weak-coupling regime of cavity QED, where damping parameters such as the pure 
dephasing rate, the spontaneous emission rate, and the cavity decay rate overwhelm the coherent coupling between the cavity 
H and the qubit. 

Here, we report that high-fidelity entanglement can be generated between homogeneously broadened matter qubits even 
in the weak-coupling regime of cavity QED. Remarkably, both spontaneous emission of the qubit and detuning between the 
photon and the qubit suppress environmental noise even for low-Q cavities, and enable distributed quantum computation to be 
performed using much less demanding technology than previously. Moreover, since appropriate use of detuning has the potential 
to overcome the huge homogeneous broadening that is the main obstacle in using NV centers at high temperatures, our analysis 
provides the possibility of using NV centers for EOs at much higher temperatures than those of current experiments [fnl l22tl . 

An outline of the proposed scheme is as follows. The matter qubit is the two ground states (|0) and |1)) of an L-type three-level 
system confined in a two-sided cavity. |0) is optically inactive, whereas |1) is radiatively coupled to an excited state |e) that is 
subject to level fluctuations due to environmental noise. Two such qubits in cavities are placed symmetrically in a Mach-Zehnder 
interferometer (Fig. [T). Initially, both qubits are prepared in ( |0) + 1 1) ) /\/2 and a single photon tuned to the cavity frequency is 
input from the left port of the first beam splitter (BS1). The state vector of the system is given by 

4(|00) + |01) + |10) + |ll))/2, (1) 

where \mn) = |m)i|n)/? denotes the two-qubit state vector and a} L (a R ) creates a photon in the left (right) path. The beam 
splitters divide a photon as a\ (a R + ia L )/\/2 and a R — > (a\ + ia R )/y/2. For the interaction between the photon and 
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FIG. 1 : Schematic view of the optical circuit. A single photon is split by a beam splitter (BS 1 ) and is sent to cavities that confine matter qubits 
with an L-type structure. After interacting with the matter qubits, the photon is combined by another beam splitter (BS2). When the photon 
reaches the target port and the detector clicks, entanglement is generated between the remote matter qubits. 

the qubit, when the qubit is in |0) (empty cavity), the input photon is perfectly transmitted through the cavity due to resonance 
tunneling. In contrast, when the qubit is in |1), the matter qubit modifies the transmitted photon. For example, as we show later, 
the matter qubit may completely prevent transmission of the photon under some conditions. Then, the photon-qubit interaction 
removes the terms ajjlO), ajjll), a^|01) and ajj|ll). In other words, the qubit state |1) acts as a "bomb" in the interaction-free 
measurement l23ll in this case. After the photon passes through the second beam splitter (BS2), the state vector is given by 

- 4=4|&> + ^4|00> + 4=ak|&>, (2) 

where \4> t ) = (|01) - |10))/\/2 and \<j>\) = (|01) + \10})/V2. A photodetector is set to count the photons that exit from the 
left port of BS2. The detector clicks with a maximum success probability of 1 /8, and the two qubits are then projected onto the 
target entangled state, \<pt). 
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FIG. 2: Schematic view of the cavity QED system that we adopted. It consists of a matter qubit and a cavity. The incoming (outgoing) photon 
fields are denoted by a, n and a' in (a out and a out ). Three types of damping are considered: environmental pure dephasing (7 P ), spontaneous 
emission (7), and cavity photon leakage (k). 

We investigate the interaction between the photon and the qubit in a more quantitative manner. Although the master equation 
has been used in previous analyses ||20[ |2lll . it is valid in principle only when the damping parameters can be regarded as 
perturbations ll24ll . In contrast, here, we solve the Heisenberg equations of the overall system including the environment in a 
non-perturbative manner. Consequently, our results are applicable to highly dissipative cases that include the weak-coupling 
regime. We investigate a cavity QED system in which a two-level matter qubit (|1), |e)) is confined in a two-sided cavity 
(Fig. |2j. The photon dynamics for the qubit state |0) is obtained by removing the matter qubit. This system is characterized 
by the following parameters: the cavity frequency u> c , the qubit transition frequency u> q , the coherent coupling between the 
cavity and the qubit g, the cavity decay rate k, the spontaneous emission rate of the qubit to non-cavity modes 7, and the 
pure dephasing rate of the qubit 7 p . The complex frequencies of the cavity and the qubit are defined by ui c = u c — in/2 and 
uj q = u q — 1(7/2 + 7 P ). We denote the destruction operators for the cavity photon and qubit by c and er(= |l)(e|), respectively. 
Their Heisenberg equations are given by 

iuj c c - iga - % ^ n/2[a in (t) + a' in (t)], (3) 
iuj q a - igc - i^jd ln (t) 

i^T P [e\ n {t)cj + ae in (t)l (4) 
where dm and ei n denote the noise operators respectively associated with spontaneous emission and pure dephasing, and ai n 
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and a' in are the incoming photon fields toward the cavity (see Fig. |2). The outgoing field operators are given by 

a ut{t) = -a' in {t)+i^fk~j2c{t), (5) 
a' out (t) = -a m {t) + i^j2c(t). (6) 

We are interested in the transmission of a single input photon. The transmitted photon consists of elastic and inelastic compo- 
nents. Thus, the state vector evolves on transmission as 



l)^U<J\l)+Ua)J\l), (7) 



where e' denotes an environmental excitation near the qubit. As we show in Appendix Bl, the fidelity and success probability 
of our EO are maximized when the spectral width of the input photon is much narrower than the cavity linewidth (i.e., the long 
pulse limit). We thus assume the long pulse limit in the remainder of the paper. The coefficients t e and ti are then determined 
by considering the linear response to a classical continuous wave. Setting a,„ = Ee~ lUat and a' in = d, n = ei n = 0, the 
dimensionless system variables (x c = iyJn/2 (c)/aj n , x a = —Wk/2 (cr)/ai n , and x c f c = k{c^ c) /2\ai n \ 2 ) are given by 



*(7/2 + 7 P + jA) 
k(7/2 + 7 p + *A) + 2 5 2 : 
Kg 

k( 7 /2 + 7 P + iA) + 2g 
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(8) 
(9) 



[7 + 2 g 2 Re(l/Q]Rc(x c ) - gi Re(x a /0 
k7 + 2 5 2 (k + 7)Rc(1/0 

where A = w 9 — w c is the detuning between the qubit and the input photon and £ = (n + j)/2 + j p + iA. t e and ti are related 
to the amplitude and flux transmissivities by t e = (a out )/cii n and \t e \ 2 + \ti\ 2 = {a\ ut a out ) /\ai n \ 2 . Thus, we have 

t e = Xc, (11) 

\U\ 2 = x c t c -\x c \ 2 . (12) 

We can confirm that inelastic transmission originates from pure dephasing, since ij = when 7 P = 0. The photon dynamics for 
the qubit state |0) is obtained by taking the g — >• limit, where we can confirm that t e = 1 and ti = 0. Therefore, the counterpart 
of Eq. © is 

a + |0) a + |0). (13) 

Using these rigorous photon-qubit interactions [Eqs. (|7) and (I13H . we reconsider the time evolution of the initial state vector 
[Eq. Q}]. Since environmental excitation inhibits photon interference at BS2, the state vector that clicks the detector is given by 

IV) = 4 (^p W + 7f «^ 2 > - 7H l<fe) ) ' (14) 

where \4> 2 ) = (|01) + |ll))/v / 2and |^ 3 ) = (| 10) + |ll))/\/2. The click probability V, reduced density matrix p, and fidelity 
F are respectively defined by V — (ip\ip), p= Tr a ,b{\'4') (V 1 !} / '(^1^) > an( l = (^t \p\4>t)- J~ and are given by 

|l-f e | 2 + N 2 /2 

|l-t e | 2 + 2|t ! | 2 ' U0J 
V = |l-t e | a /8+|*i| 2 /4. (16) 

We first examine the effects of homogeneous broadening by assuming that both detuning and spontaneous emission are absent 
(A = 7 = 0). In this case, the transmission probability through the cavity is given by |t e | 2 + |ti| 2 = kj p / (Kj p + 2g 2 ). Therefore, 
when Kjp/g 2 <C 1, the cavity nearly completely suppresses transmission of the photon and the present scheme functions with 
a high fidelity. To achieve T > 0.9 (0.95), K-f p /g 2 should be less than 0.15 (0.07). Consequently, high-Q cavities satisfying 
n lvl 9 1 ^ 1 are required to achieve high-fidelity EOs under a large homogeneous broadening. This is qualitatively consistent 
with another scheme that employs resonant input photons lfl9ll . 

Spontaneous emission usually degrades the figure of merits of quantum devices. In contrast, spontaneous emission makes our 
protocol more robust against environmental noise and relaxes the cavity conditions, so that a high-fidelity EO becomes possible 
between homogeneously broadened matter qubits even in the weak-coupling regime (g < k, 7, 7 p), as we show in Appendix B2. 
The origin of infidelity here is inelastic scattering (i.e., entanglement with the environment) that occurs while the matter qubit is 
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being excited. Spontaneous emission reduces the lifetime of the excited state and thus hinders inelastic scattering. However, in 
actual experiments, it is difficult to artificially increase the spontaneous emission rate and thus this does not provide a practical 
solution. Therefore, we look for another way to suppress environmental noise using existing technology. We consider the use of 
detuning. 

We briefly explain the physical mechanism of an EO employing a detuned photon. When there is large detuning A, Eqs. (O 
and (fTTT > give t e ~ e~ l6 , where 9 = g 2 /An. Namely, when the qubit state is |1), the input photon acquires a phase shift that 
is determined by the product of the dispersive interaction (g 2 /A) and the cavity photon lifetime (k^ 1 ) ll20ll . This mechanism 
contrasts with that of resonant cases (A = 0), where the transmitted wave is attenuated (t e < 1) through scattering or reflection. 
The fidelity can be drastically improved by detuning A because detuning hinders real excitation of the matter qubit and the 
resultant inelastic scattering. Figure [3] shows a plot of the fidelity (J 7 ) and the success probability (P) as functions of k and A, 
assuming 7 = 7 P = 2g. The cavity condition for achieving T = 0.9 is k = 0.59g when A = 0. However, this condition is 
relaxed to k — 2g by setting A = 9g. Surprisingly, high-fidelity entanglement generation is possible between homogeneously 
broadened matter qubits even in the weak-coupling regime satisfying g < k, 7, 7 p . Figure |3jb) shows that detuning reduces 
the success probability. Namely, there is a trade-off between the fidelity and the success probability. However, the success 
probability is V = 0.13% when n = 2g and A = 9g, which is sufficiently large for practical use. The dark count rate is typically 
less than 10 -7 per nanosecond so that this success probability can exceed the dark count rate even within current technology. 




FIG. 3: Contour plots of (a) fidelity and (b) success probability as functions of the cavity decay rate k and the detuning A in units of the cavity 
coupling strength g. 7 P = 7 = 2g. Even when k = 2g, high-fidelity operation (J~ — 0.9) is possible by setting A = 9g. The success 
probability will then be V = 0.13%. 



Finally, we describe a possible experimental realization of our scheme using NV centers. In a cavity QED setup composed of 
a diamond NV center and a microtoroidal cavity, the p arameters g, k, and 7 have comparable values of the order of tens of MHz 
l25ll . while 7 P is highly sensitive to temperature ill 8ll . The linewidth of the NV center will be almost lifetime limited and thus 
7 P will be negligible at low temperatures such as 7 K, whereas -f p will dominate the other parameters at higher temperatures. 
At a low temperature Qy p = O.lg and 7 = g), T = 0.96 can be attained with V = 0.13% even by a low Q cavity (k = Ag 
and A = 5g). Thus, using our scheme, it should be possible to realize an EO using current technology. Moreover, even at 
higher temperatures, it should be possible to perform an EO by our scheme with modest requirements that are expected to be 
achievable in the near future. Here, we set the parameters as 7 P /27r = 300 MHz (which corresponds to a temperature of about 
30 K 0J]), 7 /2tt = 20 MHz, g/2n = 250 MHz, k/2tt = 150 MHz, and A/2tt = 3 GHz. Entanglement can be generated 
with T = 0.90 and V = 0.96%. In principle, once this amount of remote entanglement is achieved between distant nodes, 
one can realize scalable distributed quantum computation by using purification techniques inside the nodes l26ll27ll . Therefore, 
distributed quantum computation may be possible at temperatures of tens of kelvins, which can easily be generated without 
using liquid helium I28, [29ll . 

In conclusion, we performed a non-perturbative analysis of an EO using a detuned photon as a mediator between optically 
active matter qubits. We demonstrated that this scheme is extremely robust against environmental noise so that entanglement can 
be generated between homogeneously broadened matter qubits even in the weak-coupling regime, where damping parameters 
overwhelm the coherent coupling between the cavity and the qubit. This result is particularly relevant for realizing distributed 
quantum computation by using NV centers at high temperatures of the order of tens of kelvins. Our scheme provides a practical 
way to overcome the main obstacle of using NV centers at high temperatures, namely large homogeneous broadening. 

The authors thank H. Kosaka and W. J. Munro for helpful discussions. This work was supported in part by the Funding 
Program for World-Leading Innovative R&D on Science and Technology (FIRST), KAKENHI (22241025, 23104710, and 
22244035), SCOPE (11 1507004), and NICT Commissioned Research. 
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Appendix A: cavity-QED analysis of single-photon dynamics 
1. Hamiltonian and initial state vector 

We present here mathematical details on time evolution of a single input photon in the proposed optical circuit. To begin with, 
we analyze transmission of a photon through a cavity. The physical setup is illustrated in Fig. [4] It is composed of (i) a matter 
qubit, which has three levels (|0), |1), |e)), (ii) a two-sided cavity, (iii) leak fields from the cavity (b and b' fields), (iv) noncavity 
radiation modes (d field), and (v) environmental modes causing pure dephasing of the qubit (e field). Since the state |0) is 
optically inactive, we may regard the qubit as a two-level system (|1), \e)) when investigating its optical response. Putting 
h = c = 1, the Hamiltonian is given by 



+ / dk kb\b k + y/n/Ait^bk + b\c) 



dk \kb' k % + y/ K /4ir{Jb' k + b' k f c) 



+ dk kd\d k + v /7/27r(fT t d fc + d\a) + dk ke\e k + y^o^/W 1 a{e k + ej.) 



(Al) 



where er(= |l)(e|) and c are the destruction operators of qubit and cavity photon, and a k (a = b, 6', d, e) is the destruction 
operator of a field in the wavenumber representation. The meanings of the parameters are given in the main text. The field 
operator in the real-space representation is defined by a r = (2tt)^ 1 ^ 2 J dke lkr a k . The r < (r > 0) region corresponds to the 
incoming (outgoing) field. 

At the initial moment (t = 0), we assume that a single photon is input from the b field and all other components are in their 
ground state. The initial state vector is then written as 



= / drf(r)bl\l), 

where /(r) is the wavefunction of the input photon. It is assumed to be 

f(r) = y/2/l 0(-r) exp(^w p r + r/l), 
where 6(r) is the Heavyside step function. Namely, the input photon has a pulse length I and a central frequency u> p . 



(A2) 



(A3) 




FIG. 4: The cavity-QED setup considered. A matter qubit is confined in a two-sided cavity, and a single photon is input from the left-hand 
side. 



2. Heisenberg equations 



From the Hamiltonian of Eq. (IAU . the raw Heisenberg equation for b k is given by db k jdt = —ikb k — i ^k/4tt c. This can be 
formally solved as b k (t) = b k (0)e~ lkt — i\J k/Att J* dTc(r)e~' lk ^ t ^ T ^ . As the Fourier transform of this equation, b r (t) is given 
by 

b r (t) =6 r _j(0) - i^f^j26{r)6{t - r)c{t - r). (A4) 
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Similarly, we have 

b' r (t) = b' r _ t (0) - i^/29(r)9(t - r)c(t - r), 

d r (t) = d r - t {<S)-i^e{r)0{t-r)a{t-r), 

e r (t) = e r _ t (0) - 1^27^6* (r)6»(t - r)cr t (t - r)a(t - r). 

These equations are known as the input-output relations. The Heisenberg equations for a and c are given by 

d 

—r& - 

dt 

d 



dt 



iZ) q o - ig(l - 2a^a)c - i*fy(l - 2cr t cr)<i_ t (0) - iy/2^\eL t (Q)a + ere_ t (0)], 
iuicc - iga - i y^/2\b- t (0) + b'- t (0)}, 



(A5) 
(A6) 
(A7) 



(A8) 
(A9) 



where Lo q = co, 



1(7/2 + 7 P ) and lo c = uj c — in/2. 
In the main text, the input and output fields {ai n , a' in , a out , a' out ) are defined as shown in Fig. 2. They are related to the b 
and b' fields as ai n {t) = 6_ t (0), a' in (t) = b'_ t (0), a out {t) = b +n (t) and a' out (t) = b' +0 (t). After making these replacements, 
Eqs. (3)-(6) of the main text are derived. 



3. Correlation functions 

We discuss here the following one-time correlation functions, a q {t) = {l\a{t)cr\l), a c (t) = (l|c(t)crt 1 1), (3 q (t) = 
(l\a(t)\$ in ) and (3 c (t) = (l|c(i)|* ira ). Their initial conditions are given by a q (0) = 1 and a c (0) = /3 9 (0) = /3 C (0) = 0. 
From Eqs. d A81 > and JA91 >. their equations of motion are given by 



d 


a q (t) 






-ig 




a q (t) 


di 


u c {t) 




-ig 


-iuj c 




a c (t) 


d 


' PS)' 




—iui q 


-ig 




'PS) 


di 


_ PS) _ 




-ig 


—iuj c 




. PS) 



(A10) 
(All) 



We denote the Laplace transform of a q (t) by C aq [z) = J" °° dte zt a q (t). Then, the Laplace transforms of the above equations 
are given by 





-*vT/(-*) 



C aq {z) 

Ca c (z) 
^S) 

US) 



1 



z + iuj c 

-ig 



(z-X 1 )(z-X 2 ) 
—i\J k/1 
(z-Ai)(z-A 2 )(z-A 3 ) 
2 



-i>9_ 
z + iuj q 



(A12) 
(A13) 

l/l — iu>p. The one-time correlation functions 

Next, we proceed to discuss the two-time functions such as /3 q 2 \ti,t 2 ) = (l\<r(ti)a^ (t 2 )a(t 2 )\^ in) and (3^ (t\, t 2 ) = 
(l|c(ii)<7^ {t 2 )a(t 2 )\^ in), where t\ > t 2 . Their equations of motion with respect to t\ are the same as Eq. ( lAlOt . and the initial 

(2) (2) 
conditions (ti — > t2) are given by P q '(t2,t 2 ) = PS2) and fc'(t 2 ,t 2 ) = 0. Therefore, we have 



where Ai and A2 are the two roots of (z + iuj q )(z + iuj c ) + g 2 = and A3 
are obtained by analyzing the poles of the above Laplace transforms 



P q 2) (h,t 2 ) = a q (t 1 -t 2 )PS*)i 
P (2) {h,t 2 ) = a c {tx-t 2 )P q {t 2 ). 

Repeating the same logic, general multi-time functions are written as the products of one-time functions as 

P q n) {tx,--- ,t n ) = a q (t 1 -t 2 )a q (t 2 -t 3 )---l3 q (t n - 1 -t n ), 
Pi n \t u - ..,tn) = a c (t 1 -t 2 )aS 2 -t 3 )---l3 q (t n - 1 -t n ). 



(A14) 
(A15) 

(A16) 
(A17) 



4. Wavefunctions of transmitted photon 

After interaction with the qubit-cavity system, the input photon is reflected into the b field, transmitted into the c field, or 



scattered into the d field. Time evolution of the input photon is determined by |^(i)) 



-iUt 



The state vector of the 
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transmitted component of photon is written as 



l*c(*)) 



drg (r,t)b'J + / drdx 1 g 1 (r,x 1} t)b'Je Xl 



ID 



(A18) 



Note that < r < x\ < ■ ■ ■ < t. go describes the elastic component, whereas g n (n > 1) describes the inelastic component that 
is entangled with the environmental modes (ej,). We can determine go as follows: go(r,t) — (l\b' r \^f c (t)) = (l|6j.(t)|^i n ) = 
—iy/K/2/3 c (t — r), where Eq. (IA5t is used to derive the last equality. Repeating the same arguments, we have 

g (r,t) = -iyfi/20 e (t-r), (A19) 

ffi(r,xi,f) = Hv^X-iv^/^-ziKfci - r), (A20) 

g n (r,X!, ■ ■ ■ ,x n ,t) = (-i^/Kj2)(-i^ffip) n j3 q (t - x n )a q {x n - x n -i) • • • a c (xi - r). (A21) 

On the other hand, when the qubit is in |0), the photon does not interact with the qubit and inelastic processes are absent 
accordingly. The state vector of the transmitted photon is then written as 

drg (r,tjb?\0), 



l*c(*)) 



(A22) 



where g (r, t) = lim s ^ 9o(r, t). 



5. Fidelity and success probability 

Here we investigate the density matrix of matter qubits after an entanglement operation. Throughout this section, we 
denote the photon field operator in the left (right) arm of the interferometer by aj> (ctR r ). The initial state vector is 
|^) = 2- 1 J drf(r)a{ r [\00) + |01) + 1 10) + |11)]. The beamsplitters divide a photon as a[ r -> (ia{ r + a Rr )/V2 and 
a Rr ~ * (*°Lr + a Lr)/v / 2- an d the qubit-cavity system transforms a photon as Eqs. (I A 1 8t — (I A221 > . When the photon is output in 
the left port of BS2, the state vector of the overall system is given by 

= -±=Jdr[g (r,t)-g {r,t)}a{ r \<f> t ), (A23) 
1 °° f 

- —= ^2 / drdxi ■ ■ -dx n g n (r,xi, ■ ■ ■ ,x n ,t)a f Lr e Rxi ■ ■ ■ e Rx J<j> e i), (A24) 

* n— 1 
1 °° f 

+ ~7fX] / drdx i ' ' -dx n g n (r,xi, ■ ■ ■ ,x n ,t)a) Lr e\ xi ■ ■ ■ e[ x j0 e2 ), (A25) 

v n— 1 

where \<f> t ) = (|01) - 1 10) )/y/2 is the target entangled state, \<j) el ) = (|01) + |11))/V2, and |0 e2 ) = (| 10) + \ll))/y/2. 

The success probability V of the entanglement operation, namely, the probability to click the detector, is given by V = 
(V'ilV'i)- Denoting the norm of a function / by Af(f), we have 

oo 

V = Af(g - go)/8 + X)JV(ff„)/4. (A26) 



The reduced density matrix p of matter qubits is defined by p = Tr ai e|V'i)(V'i|/ (V'ilV'i)- Therefore, 

_ M(!J0 - g0)|&)(&| +Sn=l^(fln)(|^l)(0el| + |^e2><0e 2 |) 

P ^(5o-fl«>) + 2E^i^(fln) 

The fidelity J 7 between p and the target state |<^ t ) (<^> t | is given by 



(A27) 



AA(ff -.9o) + Er=i2AA( ff „) ■ 1 8j 

It is of note that the infinite sum of N(g n ) can be carried out analytically. Using the Laplace transforms of |a c | 2 , \a q \ 2 

and |/3 9 | 2 , we have 

In the long pulse limit of I — s- oo, Af(g — go) and E^Li A/"(<? n ) respectively reduce to |1 — t e \ 2 and l^l 2 as discussed in the 
main text. 
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Appendix B: numerical results 

In this section we present the numerical results that are not presented in the main text. We assume oj p = oj c throughout this 
section and denote the qubit-cavity detuning io q — uj c by A. 

1. Pulse Length 

First, we observe the effects of a finite pulse length I of an input photon. Assuming a dissipation-free (7 = j p = 0) and 
resonant (A = 0) case, the success probability V is plotted as a function of I for several values of k in Fig. [2 a). We can observe 
there that V becomes independent of I for I ^ k" 1 and reaches the limit value given by Eq. (16) of main text. This implies 
that the long-pulse limit, where the input photon can enter the cavity perfectly, is achieved when the spectral width of input 
photon is much narrower than that of cavity (k). For shorter pulses, the cavity filters out the off-resonant components of 
input photon and the success probability is decreased accordingly. In the short-pulse region, the success probability becomes 
proportional to I since it is determined by the overlap between the spectra of input photon and cavity. 

Figure [2b ) shows the ^-dependence of fidelity. As expected, T becomes independent of I in the long-pulse limit and the 
limit value is given by Eq. (15) of the main text. However, in contrast with Fig. [2a), the fidelity is insensitive to I also in the 
short pulse region. This can be understood intuitively as follows. Once the photon enters the cavity, its property is determined 
by the cavity linewidth and becomes irrelevant to the original linewidth determined by I. We can observe that both the success 
probability and the fidelity are maximized in the long pulse limit. 



K/g=0.2 



K/g=2 



K/g=2cr 

io~ 2 to" 1 1 10 10 2 
/ (units of 1/g) 

FIG. 5: Dependences of (a) success probability and (b) fidelity on the input pulse length I. 7 = 7 P = and A = 0. The values of k are 
indicated in the figure. 



2. Spontaneous Emission 

Here we observe the effects of nonzero 7. Assuming a noisy environment Qy p = 2g) and a resonant input photon (A = 0), 
the fidelity T is plotted as a function of k and 7 in Fig.[6|a). We can confirm that the cavity condition is substantially relaxed 
by a nonzero 7. In order to achieve T = 0.9 for example, k = 0.08g is required when 7 is absent, whereas this condition 
is relaxed to n = 0.59<? when 7 = 2g. Usually, spontaneous emission into irrelevant modes leads to dissipation of quantum 
devices and lowers their figure of merits. However, this is not the case with the present scheme. The origin of infidelity here is 
inelastic scattering (in other words, entanglement with environment), which occurs while the qubit is being excited. Spontaneous 
emission makes the lifetime of excited state shorter and thus hinders inelastic scattering. The success probability V is shown 
in Fig.[6|b). It is observed that V is lowered by 7. Thus, a high-fidelity operation becomes possible at the expense of a lower 
success probability. 
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